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I. INTRODUCTION 


In this paper we discuss important relativistic effects and issues which must 
be considered in the interpretation of current measurements such as ranging 
measurements to LAGEOS and to the moon, in the implementation of the Global 
Positioning System, in the synchronization of clocks near the earth's surface, and in 
the adoption of appropriate scales of time and length for the communication of 
scientific results. 


II. FEATURES OF GENERAL RELATIVITY 

We restrict our discussion to Einstein's General Theory of Relativity (GR). GR 
describes events in a four-dimensional, space-time manifold (with coordinates X^), in 
which the invariant interval ds between events separated by dX^ 1 is given in terms of 
the metric tensor G^v by: 


- ds 2 = GjxvdX^ dX v (summed on repeated indices). (1) 

The interval ds has the following interpretations. First, a test particle in free 

fall, such as an earth-orbiting satellite, follows a geodesic path, along which ds is an 

extremum. Thus, equations of motion of bodies in free fall can be derived from a 

knowledge of the metric tensor. Second, the proper time elapsed on a standard clock 

(such as an atomic clock) in free fall through the interval dX* 1 will be given by I ds I . 

Third, a pulse of electromagnetic radiation will travel along a null geodesic, ds = 0; 

this is another way of stating the constancy of the speed of light. 

Also, ds is a scalar quantity, invariant with respect to arbitrary 

transformations of coordinates. Much of the calculation entailed in working out the 

implications of GR involves finding coordinates in which observations can be readily 
interpreted. Thus, along with the knowledge of a metric tensor, in a particular 
coordinate system, goes a procedure for interpretation of the theory. 

III. BARYCENTRIC COORDINATES 

Except for small structure effects to be mentioned later, solar system bodies 
can be described by an approximate point mass metric derived from Einstein's field 
equations by Eddington and Clark (1938). The metric tensor can be written in the 

following interesting form, with Goo accurate to order (V 4 ): 
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where the retarded time is T — T — | X X^(T') | /c. Thus, in the radiation gauge used 

for the solution of the field equations by Eddington and Clark, a retarded Lienard- 
Wiechert potential (Jackson 1975) appears. The other components are: 
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where X A , V A , and A a are the position, velocity, and acceleration of the A 1 * 1 mass, 
and a prime on a summation symbol means that terms which are indefinitely large 
are to be omitted. M A is the Schwarzschild mass parameter of the A th mass; velocities 

are measured in units of c. We make use of the abbreviation R = | X A " X B I which 

is a function of X because both X A and Xg depend on X . We introduce the 
"normalized" or external part of the metric by defining the negative of the potential 
in the neighborhood of the earth's mass Mg due to external sources (the subscript E 
denotes quantities describing the earth): 


u (e) =r 


A |X- X A 


"= X' M A /R A ; U e =U (e) (x=X E ) . 


Using retarded time, expanding the retarded potential gives three retardation 
corrections in addition to the static potential: 
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where the unit vector N A is given by: N =(x _ X ) / | X"X | 
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One can see the origin of terms of interest arising from retardation of gravitational 
signals in this coordinate system. Observation of such terms would be of great 
interest as it would provide indirect evidence for the existence of gravitational 
waves. 


The equations of motion have been derived in detail by Moyer (1981 a,b) and 
others and may be written: 


k / 


M.X n . 

v « A BA 

B t 




1'4U .+V^ - - 4V r • V 

c 


a t, a - R ce "’ b b a 


|( X AB -Vj 


R 


AB 


1 

2 X AB * A B 


+41' 

B 


m b( x ab • v A )( 


V* - v k B ) 


R 


AB 


( 6 ) 


X AB- V b)( V A V b) 7 
+ 2 

B K AB B AB 


where: 


X AB - X A - X B • U A - u<e> ( x = x a) • Nab = X AB I R AB ■ 


(7) 


For the case of an earth-orbiting satellite such as LAGEOS, significant perturbations 
arise from the term with the coefficient of -3/2 (which arises from retardation in 
this picture). This particular term could give rise to accelerations as large as 4.1 x 
10 -8 meters/sec 2 . The Newtonian solar tidal accelerations of LAGEOS are only 10 5 

larger. 


IV. LOCAL INERTIAL FRAME 

It is natural to do numerical computations of planetary ephemerides in solar 
system barycentric coordinates, using the above equations of motion that include 
effects due to all solar system bodies and to post-Newtonian relativistic effects. 

However, most observations are made from the earth — a locally inertial, freely 
falling platform — and do not reveal the existence of such effects. The proper 

interpretation of barycentric coordinates in terms of measurable quantities is key. 

Our approach is to construct a transformation of coordinates from barycentric 
to local inertial coordinates — normal Fermi coordinates (Manasse and Misner 1963; 
Ashby and Bertotti 1984, 1986). In these coordinates, apart from the gravitational 

effects due to the Earth's mass itself, one can demand that time be measured by a 
standard clock near the position where the observations are made. In contrast, the 

coordinate time variable in the Eddington-Clark metric is measured by a standard 
clock at rest at infinity. Also, one can demand that the spatial coordinates be 

measured by standard rods or, in other words, that lengths and times be related in 
such a way the electromagnetic signals propagate with the defined value of the speed 
of light, c = 299, 792, 458 meters/sec. 
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Referring to Figure 1, X fA (s) specifies a world line G, a timelike solution of the 
geodesic equations of the mass Mg in the metric given by Eqs. (2 and 3), with infinite 

self-interaction terms omitted. Four mutually orthonormal vectors A ( J ,a= 0,1, 2,3 are 

introduced which are carried parallel to themselves along G. The zero lh member of 

the orthonormal tetrad is the tangent vector to the geodesic. Given a field point P(X H ) 
near G, a spacelike geodesic S is constructed which passes through P and intersects G 
orthogonally at proper time = s as measured on a standard clock falling along G. 
The tetrad forms the basis for coordinates in the (almost) locally inertial frame. The 
first stage of the construction is finding these basic vectors. 

The second stage is calculation of the transformation equations themselves by 
means of: 


X^(P) = (P 0 ) + A J 


2 


(P ) A “ A* 5 x 1 x J 
v 0' (i) (j) * x 
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X* x j x k + 0(x 4 ) 


( 8 ) 


Lower case letters x 1 are used to denote coordinates in the local frame. The 
coefficients of x 1 in Eq. (8) are evaluated at P Q and are functions of s (or x°) only. The 

coefficients of the third and higher order terms in Eq. (8) are obtained from the 
equation of the space-like geodesic S. Cubic and quartic terms are necessary in order 
to verify that the field equations are satisfied in the local frame. 

The leading terms in the resulting coordinate transformations have simple 
physical interpretations and can be written as follows: 


X° = |K dx° + v E -7 = X° (P Q ) + v E -7 +... 

X k = X^(X° <p 0 » ♦ x‘( 1 - u e - A e -7 - \ U e mn x™ x") 

+ 2 V E 'VO + !l kj S jm x m + j ? A k +... 

In the first term in Eq. (9), the factor K is given by: 


K -1 = ds/dX° 


G = V(-O iiv dX'‘dX”) 


u .- v e/2 


(9) 

( 10 ) 


( 11 ) 


and gives the rate at which proper time elapses on a standard clock falling along 
with the origin of coordinates, with respect to barycentric coordinate time. Both 
special and general relativity effects are incorporated in K. U e , the negative of the 
potential at the position of the falling clock due to all sources except the earth, gives 
rise to a gravitational frequency shift, and the term in Vg represents the second- 
order Doppler shift of the moving clock. 
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The second term in Eq. (9) arises from the well-known breakdown of 
simultaneity. This term plays an extremely important role in the coordinate 
transformation. The factor (1-U e ) in Eq. (10) represents an overall length scale 

change, due to the external gravitational potential. The term quadratic in the Earth's 
velocity in Eq. (10) represents Lorentz contraction, and the term involving the 

antisymmetric quantity represents a slow rotation (geodetic precession) of the 

inertial frame axes with respect to distant stars (of about 19 milliarcsecs per year). 


To describe gravitational effects in the local frame, one must perform a tensor 
transformation of the components of the metric, where the necessary partial 

derivatives, 6X ct /8x Ji , are obtained from the coordinate transformations. One must also 
substitute the transformations into the expressions upon which G^ v depends. For 

i —* 0 I . ... 

example, in the quantity X X^ (X ) , which occurs in the earth’s potential, it is 

—> —■ * o 

assumed the field point X and the source point Xg (X ) are evaluated simultaneously 

in barycentric coordinates. The substitution must be done with exceptional care 
because of the relativity of simultaneity. In performing this substitution, one finds 
that Lorentz contractions and the relativity of simultaneity give rise to similar 
corrections, but that due to simultaneity is twice as large at that due to Lorentz 
contraction and is of the opposite sign. One finds: 


2M r 


x- X E (X°) 


^{, + U c «.V2 + (7.V E / r ) 2 /2} . 


( 12 ) 


The term in U e arises from rescaling of lengths, the term in V £ from a combination of 
simultaneity breakdown and Lorentz contraction, and the term in A p from spatial 
curvature (quadratic terms in the transformation of coordinates). 

There are a number of additional contributions to g 0Q which are proportional 
to 2M E /r. These are as follows: 

(1) An overall multiplicative factor contributes an additional 

(K 2 - 1) x 2M E /r - (2U + V 2 ) 2M E /r; 

(2) Nonlinear velocity correction terms to the masses which remain in G 00 

2 

in addition to the velocity terms in the retarded potential of 3M E V E /r; 

(3) A cross-term between the external potential and the earth's potential in 
the squared potential term in G 0Q of -4M E U e /r; 

(4) Nonlinear interaction terms in G 00 of -2M E U e /r; 

(5) Contributions from "magnetic" terms, Goi of -8M E V E /r; and 

(6) Contributions from G . of +2M C V 2 / r . 

ij E E 

Collecting all the terms, we obtain: 
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(13) 
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This result must now be combined with the expansion of the retarded potential, Eq. 
(5). It is easily seen that all corrections cancel, leaving just the static potential term 

2Mg/r in the local frame. Thus, what appears to be a retarded potential from a 

moving source in barycentric coordinates appears as a static potential in a local 
inertial frame falling along with the mass. 


Obtaining the metric tensor in the local frame by transformation is 
complicated; it has been illustrated above for a term of one type. This has only been 
carried out for a restricted class of models because of the complexity of the 
calculations. In the general case, one can calculate the terms linear in x k in the local 

frame and Find that there are about 3 dozen terms, of 16 different types, all of which 
cancel out. This is as one would expect on the basis of the principle of equivalence, 
according to which gravitational forces due to distant bodies can be transformed 
away locally (i.e., at the origin of local coordinates), by transforming to a freely 
falling inertial frame. 


For a model in which one considers the earth to be falling around the sun (of 
mass Mg) in a circular orbit, the terms quadratic in local spatial coordinates can be 

calculated. The contributions to g 00 consist of: a Minkowski term, the Newtonian 

potential term, a contribution from the nonlinear Schwarzschild field of the earth, 

Newtonian solar tidal terms, nonlinear earth-sun interaction terms, nonlinear solar 

tides due to interaction between the sun and the earth, and higher order solar tides. 
The spatial part of the metric in the local frame has the usual spatial curvature term 
due to the earth, plus solar tidal corrections. The g Q . terms contribute magnetic 

effects. 


Calculations have also been performed using the PPN metric with parameters 
P»Y*£1> and £2 (Shahid-Saless and Ashby, in preparation). These give rise to some 
interesting effects involving shifts of the center of mass of the earth-sun system. 
The elliptical orbit case of the earth has also been treated, assuming, however, that 
the earth can be treated as a test particle so certain types of nonlinear interactions 
can be neglected. To verify the cancellations in more realistic cases requires 
considerable calculation. Computer algebra programs help some, but there is a 
tendency for such programs to fill up memory and fail in working on this problem. 

V. EQUATIONS OF MOTION 

Equations of motion can be derived from the local metric in a straightforward 
way. Here, we shall give only estimates of the orders of magnitude of the resulting 
accelerations of an earth-orbiting satellite caused by post-Newtonian effects. In 
evaluating the orders of magnitude, it is assumed that the satellite orbit satisfies 
v 2 - c 2 M E /r. For LAGEOS r/10 9 cm - 1.2. 
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TABLE 1 


ESTIMATES OF ACCELERATIONS OF AN EARTH SATELLITE DUE TO RELATIVITY THE MAGNITUDES OF THE 
ACCELERATIONS DUE TO NEWTONIAN FORCES OF ATTRACTION TO EARTH, AND THE SOLAR TICES, ARE GIVEN 
FOR COMPARISON (N. ASHBY AND B. BERTOTTI 1984). 


SOURCE 


MAGNITUDE 

Newtonian potential 

c 2 M E /r 2 

— 4x10 2 [(10 9 cm)/r] 2 cm/sec 2 

Solar tides 

c 2 M(y/R 3 

— 4xl0~ 5 [r/(10 9 cm)] cm/sec 2 

Nonlinear earth field 

c 2 M 2 /r 3 

E 

— 2x 10" 7 [(10 9 cm)/r] 3 cm/sec 2 

Nonlinear solar tides 

c 2 M 0 2 r/R 4 

— 4xl0" 13 [r/(10 9 cm)] cm/sec 2 

Solar ’’magnetic terms” 

c 2 7[MQ 3 M E r/R 7 ] 

— 8xl0~ 13 V[r/(10 9 cm)] cm/sec 2 

Earth-sun interaction 

^MqMe/R 3 

— 2x10" 14 cm/sec 2 


Note the relativistic perturbations are many orders of magnitude smaller than first 
calculated in the barycentric frame. This is because of the cancellations which 
occur upon transforming to a freely falling coordinate system. Rubincam (1977) has 

studied the effect of the nonlinear earth field. 

The above relativistic orbital perturbations are so small that they cannot be 
expected to significantly affect LAGEOS’s orbit. However, the nonlinear solar tidal 
term grows with distance from the earth and is much larger at the orbit of the moon. 
Combining the effects of this perturbing term on semimajor axis of the moon, 

eccentricity of the moon, and on the moon's mean motion, gives rise to a net 

perturbation on the distance between earth and moon of : 

5r — 4.4 cm x cos(2f M + 2co - 2f E ) (14) 

where f M and f E are the true anomalies of the moon and of the earth, respectively. 

This effect is unfortunately obscured by its high correlation with effects due to the 

solar tides themselves. 

A more interesting set of effects, which have not yet been fully explored 
analytically, arise from multipole contributions to the earth's field. The leading 

_ o 

quadrupole contribution is about 10 of the main monopole term and relativistic 
corrections arising when this quadrupole is viewed as moving in the barycentric 

o 

frame can be expected to be 10~ smaller, thus such relativistic accelerations can be 

estimated to be about 4 x 10~ 9 cm/sec 2 for LAGEOS, which is significant. Inclusion of 
such effects results in significant improvement in fitting the ranging data for 
LAGEOS (Ries, private communication). 

VI. CHOICE OF TIME COORDINATE— TDT 

Further transformations must be made since time standards laboratories 
having clocks used to define the SI second are on earth's surface, subject to additional 
motions due to earth rotation, and to the full gravitational potential of the earth 
including higher order multipole contributions to the potential. 
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Since the geoid is a surface of gravitational equipotential in the earth-fixed 

rotating frame, atomic clocks at rest on the geoid on the rotating earth all beat at the 
same rate. If one compares a clock near the pole with one on the equator, one finds 

that the one nearer the pole is close to the earth’s center and is therefore beating 
more slowly due to a gravitational redshift effect. However it is also closer to earth's 
rotation axis and moving more slowly due to earth’s rotation, and is subject to less 
time dilation (second-order Doppler shift). These effects cancel on the geoid. 

The rate of a standard reference clock falling along the geodesic G, with rate 
correction given by the factor K in Eq. (11), does not incorporate the effects due to 

earth’s mass and rotation. An additional correction is needed to obtain a new 

coordinate time x° SI , corresponding to the definition of the SI second. This is given 

by: 


dx 


0 


, 1 M E r . T m 1 , . ,2 

= 1 "277 (1 + J 2 /2 >-2 ( ® a l /c) 


dx 


o 


(15) 


where ai is the equatorial radius of the earth. Then, apart from periodic terms, 

.0 


dx 


SI 


M, 


dX° 


= 1 - 


2ajC 


^2 (1 + ty 2 ) ~ (“Vc) 2 # - U e - V E /2 3 1 - L 


(16) 


— 8 

where L — 1.55 x 10 ; Hellings (1986) has given a more complete discussion of this. 
Barycentric Dynamical Time (TDB) runs at the same average rate as the SI second, so 
^TDB = (1 - L)X°. Thus TDB clocks beat at the same average rate as earth-borne clocks, 

from the point of view of an observer in the barycentric system. The TDB clocks beat 
more slowly, by the factor 1-L, than coordinate clocks in the Eddington-Clark metric. 
Therefore, to maintain a universally defined numerical value for the speed of light c, 
the length unit in TDB coordinates must be physically longer than the length unit in 

EC coordinates. Then since (GM/c ) TDB represents a physical length as measured 

using a TDB meter stick, the numerical value of (GM/c 2 ) TriD will be less than it is in EC 
coordinates. But the speed of light, c, is the same in the two unit systems, so (GM) TDB 
(1-L) (GM) S] . 


An outstanding issue is the definition of TDT — "terrestrial dynamical time." 
The problem is to define a terrestrial time scale which is as closely related to TAI in 
rate as is possible. Problems include how to specify the position of the master 
reference clock for the TDT (the geocenter is preferred), how to specify the rate of 
TDT in relation to TAI, and how to initialize the TDT clock since its position is not that 
of any clock contributing to TAI. Some synchronization convention must be adopted 
to resolve ambiguities. 


VII. RELATIVISTIC EFFECTS IN THE GLOBAL POSITIONING SYSTEM 

This area of application has been discussed elsewhere (Hellings 1986) so only a 
summary of the relativistic effects is given here. For clocks in GPS satellites or on 

the earth's surface, it is useful to synchronize to agree with hypothetical clocks 
synchronized in the local inertial frame. There are three residual relativistic 
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effects: second-order Doppler frequency shifts of moving clocks, gravitational 

frequency shifts, and the Sagnac effect. 

To the readings of each atomic clock, systematic corrections can be applied 
based on the known positions and motions of the clocks, such that, at each instant, 
the coordinate time thus produced agrees with the reading of a fictitious standard 
clock with which it instantaneously coincides and which is at rest in the local 
inertial frame. GPS coordinate time is thus generated by systematically modifying 
the proper time elapsed on standard clocks. 

Standard clocks on the earth’s surface provide the reference rate. Due to 

gravitational potential and motional effects, clocks in the satellites require a 

fractional rate offset correction having the value: 

3GM P GMp - in 

f f (1 + JJ2) - (u>aJc) 2 /2 ~ - 4.465 x 1(T 10 . (17) 

2ac ajC 1 

In order for the SV clock to appear to an observer on the ground to beat at the chosen 
frequency of 10.23MHz, the SV clocks must be offset so that, to an observer in the SV 

rest frame, the frequency in 10.23 x (1-4.465 x 10 ^)MHz - 10. 22999999543MHz. 

An additional variable correction term, arising from a combination of 
gravitational frequency and second-order Doppler shifts in case the SV orbit is not 
circular, must be applied to the SV cock to yield GPS coordinate time: 

At sy = +4.4428 x 10“ 10 sin E(t) (18) 

where e is the orbit eccentricity and E is the eccentric anomaly. 

For a standard clock transported near the earth's surface, the following 
correction must be applied: 

At'=Jp ath (ds/c) [ 1-OD-V/c 2 + (v'/c) 2 /2] + 2coA’ e /c 2 (19) 

where primed quantities are measured in the earth-fixed rotating frame. This gives 
the prescription for correcting the atomic clock reading Jds/c to account 
respectively for gravitational frequency shifts, second-order Doppler shifts, and the 
Sagnac effect. For clocks near the Earth's surface, 4>-<t> 0 cl gh> where g is the 

acceleration of gravity and h is the height above the geoid. The term (v/c) / 2 
corrects for time dilation of clocks moving relative to the ground, and the last term 

expressing the Sagnac effect is equal to 2coA , g/c^, where Ag is the equatorial 

projection of the area swept out by the position vector r 1 of the clock in the rotating 
frame. 


For synchronization by means of an electromagnetic signal, ds vanishes along 
the path. The elapsed coordinate time during propagation of the signal is: 

At^jpath da' [l-(<M> 0 )/c 2 ] + 2 wA’ e /c 2 , (20) 
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( 20 ) 


At " j path d<7 ' [l-(<t>-V /c2 l + 2 ® A V c2 ’ 

where da 1 is the increment of proper distance along the path of the signal and A'e is 
the equatorial projection of the area swept out by the position vector of the electro- 
magnetic pulse in the rotating frame. In the above expressions. 


MG 

<M> 0 = J—[ 1- J 2 (a 1 /r) p 2 (cos 0)] - (cor' sin 0) 2 /2 

GM 

+ -^"*(1 + J 2 /2) + (0)aj) 2 /2 . 


( 21 ) 


In GR, the same Sagnac correction terms arise whether synchronizing clocks 
by slow transport of portable clocks or by transmission of electromagnetic signals in 
the rotating frame. 
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DISCUSSION 


NIETO: In the case that you have a monopole plus J2 representing the earth, what is 

the motion of the geocenter about the geodesic? 

ASHBY: The geocenter wobbles about a timelike geodesic with an amplitude of about a 

meter. 



FIG. 1. — Diagram showing the local inertial frame falling freely along the geodesic G. The space-like geodesic S 
is constructed by dropping a geodesic from the field point P to P 0 , that intersects G orthogonally. 
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